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Disc instability in RS Ophiuchi: a path to Type la supernovae? 
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ABSTRACT 

We study the stability of disc accretion in the recurrent nova RS Ophiuchi. We construct a 
one-dimensional time-dependent model of the binary-disc system, which includes viscous 
heating and radiative cooling and a self-consistent treatment of the binary potential. We find 
that the extended accretion disc in this system is always unstable to the thermal-viscous in- 
stabiUty, and undergoes repeated disc outbursts on ~ 10-20yr time-scales. This is similar to 
the recurrence time-scale of observed outbursts in the RS Oph system, but we show that the 
disc's accretion luminosity during outburst is insufficient to explain the observed outbursts. 
We explore a range of models, and find that in most cases the accretion rate during outbursts 
reaches or exceeds the critical accretion rate for stable nuclear burning on the white dwarf 
surface. Consequently we suggest that a surface nuclear burning triggered by disc instability 
may be responsible for the observed outbursts. This allows the white dwarf mass to grow over 
time, and we suggest that disc instability in RS Oph and similar systems may represent a path 
to Type la supernovae. 

Key words: accretion, accretion discs - instabilities - novae, cataclysmic variables - binaries: 
symbiotic - supernovae: general - white dwarfs 



1 INTRODUCTION 

The recurrent novae (RNe) are a small set of close binary sys- 
tems which are observed to have nova-like outbursts every ~ 10- 
100 years. The main reason for interest in them is their possi- 
ble relevance in the search for progenitors of Type la supernovae 
■ (SNe la). It is generally agreed that SNe la occur when a white 
' dwarf (WD) is driven over t he Chandrasekhar mass Mch (e.g., 
\ iHillebrandt & NiemeveJio og). The main obstacle to arranging this 
. is that accretion of hydrogen-rich matter on to a WD generally ap- 
pears to result in a thermonuclear runaway - a nova explosion - 
once a certain mass has been accreted. All novae recur, and RNe 
earn their name simply because their recunence times are short 
enough to be observed. However, nova explosions remove most 
or all of the mass that was accumulated since the last explosion, 
and it seems likely that the masses of accreting WDs are slowly 
eroded by this pro cess, rather than increasing (e.g.. iNomotoll 1 983 : 
lYaron et al.ll2005l) . The only exception to this is that the accreting 
mass can apparently undergo steady nuclear burning if it arrives at 
a rather narrowly confined rate Msteadv - 3 -8 x lO^'Moyr"' (e.g., 
IShen & BildsterJl200llNomoto et al.ll2007h . 

Two solutions to this problem are currently suggested. One 
is that SNe la result from the merger of two WDs with total 
mass exceeding Mch l lWebbinkllTgslllben & Tutukovlll984h : the 
other model picks out an unusual subset of accreting WD systems 
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( IWhelan&Iber]|l973l: iNomotd [1983) . Observations of the pro- 
genitors of Type la SNe are scarce, but generally point towards 
growth by accretion rat her than a double-degenerate origin (e.g., 
IVoss & Nelemansll2008h . However, it is far from obvious how this 
accretion proceeds. The vast majority of accreting WDs are cat- 
aclysmic variables (CVs), where the companion star has a lower 
mass than the WD. Mass transfer here requires orbital angular mo- 
mentum loss, or in rare cases nuclear expansion of the companion. 
However if the close companion is more massive than the WD, 
mass transfer is self-sustaining on the thermal time-scale fxH of the 
companion, a so-called supersoft source. For a likely companion 
mass ~ IMothis mass transfer rate is plausibly close to M steadv; 
allowing the WD mass to grow dLi & van den HeuvelllT997l) . The 
model requires some fine tuning to ensure that the mass transfer 
rate remains close to Mstcady. and can only increase the WD mass at 
most t o the current companion mass. R ecently Gilfanov & Bogdd^ 
1 I2OIOI) and iBogdan & Gilfanov! i I^dTH) have shown that the total 
soft X-ray emission of both early- and late-type galaxies implies 
numbers of supersoft sources which are too low to account for the 
global SNe la rate. 



The RNe offer intriguing insight into these problems. Fit- 
ting their very short recurrence times into the general picture of 
classical novae requires the RNe to have WD masses Mi very 
close to Mch- The argument runs as follows (e.g.. lTruran & Liviol 
IT98^ . A thermonuclear runaway requires a critical pressure /'(■rii ~ 
2 X 10" dyne cm"', independent of WD mass, at the base of the 
accreted layer of mass AM so that 
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GMiAM 



(1) 



where i?„d is the WD radius. The short recurrence time-scale of 
RNe allows little time to accumulate matter, so AM is unusually 
small. From Equation [T] we see that AM oc /S^^, so recurrent no- 
vae require small R^i and thus WD masses very close to Mch- 
More recent work suggests that other factors, such as the com- 
position and thermal profil e of the igniting layer, may also affect 
the critical mass AM (e.e.. IShen & Bildstenll200^ : [ Kasliwal et al.l 
but short recurrence time-scales still generally require WD 
masses close to Mch- Unfortunately the difficulty of measuring 
WD masses in CVs means that t here is little strong evidenc e ei- 
ther for or against this prediction. iThoroughgood et al.1 lioOlh find 
Ml = 1.55 ± 0.24Mo for U Sco, consistent with the theoretical 
prediction, but few other dynamical masses have been measured. 

However, other interpretations of the recurrent novae have 
been suggested. In part icular, the observation of a radio jet in the 
recurrent nova RS Oph jSokoloski. Rupen. & Mioduszewskil2008h 
is difficult to reconcile with the idea of a thermonuclear runaway 
jKing & P ringle 2009). Well-collimated jets require the presence 
of an accretion disc, whose inner parts would be destroyed by a 
near-spherical nova explosion since the nuclear energy release per 
gram considerabl y exceeds the gravitational b inding energy at the 
WD surface (e.g.. lFrank, King. & Raine"2002'). Therefore, at least 
at the onset of the outburst, thermonuclear energy release cannot 
dom inate the gravita t ional contribution. 

iKing & Pringld ^2009.) noted that the very long orbital period 
of RS Oph (~ 454 days) and the relatively low mass M-giant sec- 
ondary imply a wide orbital separation (> lO'^'cm) and a very large 
accretion disc (see also Wynn 2008 ,). The outer regions of such a 
disc are cool, and the disc is therefore subject to the usual thermal- 
viscous disc instabil ity which p owers dwarf novae (see e.g. Lasota, 
2001 for a review). I King & Pr ingle (2009) therefore attempted to 
interpret the outbursts of RS Oph as very large dwarf nova events, 
where the outburst luminosity comes solely from gravitational ac- 
cretion. The aim of this study is to test this interpretation in detail. 

In this paper we present a one-dimensional, time-dependent 
model for disc accretion in the RS Oph system. Our model includes 
viscous heating and radiative cooling, as well as a self-consistent 
treatment of the binary potential, and is able to reproduce outbursts 
due to the thermal-viscous instability. We also model the observed 
emission from the disc, during both quiescence and outburst. The 
disc is fed by Roche lobe overflow from the giant secondary, and we 
find that the disc is always unstable to the thermal-viscous insta- 
bility. The outbursts naturally recur on ~ lOyr time-scales, but the 
accretion disc is not sufficiently luminous to explain the observed 
outbursts. The total mass accreted during these outbursts can be 
large, < IO^'Mq, but much of the accretion occurs at rates large 
enough to trigger nuclear burning on the WD surface. We discuss 
the consequences of our results for the observed outbursts of RS 
Oph, and speculate as to their likely origin. 



2 MODEL 
2.1 Model 

The evolution of the gas surface density, /) in the ac- 
cretion disc around the primary is governed by the equation 
jLynden-Bell & Pringle 1974; Pringle 1981) 



'di 



RdR 



dR 



(GMi)i/2 



+ 2jnfall(/?, f) . 



(2) 



Here t is time, R is cylindrical radius, v is the kinematic viscosity 
and Ml is the mass of the primary. The first term on the right-hand 
side represents viscous evolution of the disc, while the source term 
%nfA\(R, t) represents the mass infall on to the disc. The second term 
describes the response of the disc to the tidal torque from the sec- 
ondary. A{R, R2) is the rate of speci fic angular momentum t ransfer 
from the secondary to the disc (e.g., lLin & Papaloizoull98^ . given 
by 



A(R,R2) = — 



rGM, I R 



2R 



where q = M2/M1 is the binary mass ratio and 
Ap = max(H,\R-R2\). 



(3) 



(4) 



Here R2 is the binary separa tion and and H = cJQ. is th e disk scale- 
height. We adopt a standard fshakura & SunvaevI d 1 9731) tt- viscosity 



(5) 



where cs is the local sound speed, Q = ^jGMJR^ is the Keplerian 
angular frequency, and the dimensionless parameter a specifies the 
efficiency of angular momentum tr ansport (thought to be due to th e 
magnetorotational instability, e.g.. lBalbus & Hawlevlll99lLll998l). 
We u se a simplified form for the energy equation (e.g. Cannizzol 
I1993I) 



_ . . 2(6^ - e_) 57^ 

dt cpE dR 



(6) 



Here is the midplane temperature and Q+IQ- are the instanta- 
neous heating/cooling rates. The radial velocity in the advective 
term in Equation|6]is taken to be the vertically averaged value, and 
is computed as 

3 d 



ur 



, • (7) 

The specific heat capacity Cp is approximately constant for most 
temperatures of interest, but rises dramatically at lO^'K (due to 
the ionizati on of atomic hyd rogen). We adopt the functional form 
for Cp from lCannizz 3( fl993h . The heating term 2+ is evaluated as 



2+ = Gv + Gtid + 2hs • 

Qy is the heating from viscous dissipation, given by 
9 , 



(8) 



(9) 



The heating due to the tidal torque from the secondary is evaluated 
as (e.g.. iLodato et al. 2009) 



e.id = [^2 - n(R)]A(R)l(R)dR, 



(10) 



where ^2 is the orbital frequency of the binary. The heat input at 
the "hot-spot" w here the accretion stream lands on the disc is given 
by ( lLasotall200lh 



Ghs = v 



1 



2R„, 



InARv, 



■ exp 



~AR^ 



(11) 



Here rj = 1 is an (arbitrary) efficiency factor, Minf^ii (= 
^ 2nR%nf.^\\dR) is the accretion rate from the secondary, Roui is the 
outer disc radius (computed as the radius where X falls below 1 
g cm"^), and A/?hs = 0.1/?out is the width of the heated region. In 
practice viscous heating dominates: 2hs is almost completely neg- 
ligible, and tidal heating is only significant very close the the outer 
edge of the disc (and even then is not dominant). 

We adopt a "one-zone" mod el for cooling (as used by, for ex- 
ample. ljohnson & Gammiell2003l) , so that 
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(12) 



Here, the second term allows for a smooth transition between the 
limits of optically thick and optically thin cooling, ctsb is the 
Stefan-Boltzmann constant, and the last term is prescribed in a 
manner than enforces a minimum disc temperatute Tn,;,, = lOK. 
The vertical optical depth t is computed as 



1 



(13) 



with the midplane density evaluated from the surface density by 
assuming that the disc is vertically isothermal: 



(14) 



The local sound speed is essentially a means of keeping track of 
the local pressure, and includes contributions from both gas and 
radiation pressure: 



4crsBr;* 



(15) 



Here 7? is the gas constant, ji is the mean molecular weight, and c 
is the speed of light. In practice, radiation pressure is only signifi- 
cant at high temperatures. We adopt s tandard numerical fits for both 
k{T,p) and ii(T,p): tha t of IZhu et al . (2007, 2008) for the opacity, 
and that of lHurel ( l200oh for the mean molecular weight. 



2.1.1 Numerical Method 

We solve Equations |2] and |6] using an explicit finite- 
difference method, on a grid that is equispaced in (e.g., 
IPringle. Verbunt & Wadd [l983) , using 1587 grid cells spanning 
the range [lO'^cm, lO'^cm]. We use a staggered grid, with scalar 
quantities evaluated at at cell centres and vectors (velocities) at 
cell edges. The derivative in the advective term (Equation IT} is 
evaluated as the first-order, up-winded value. We adopt zero-torque 
boundary conditions at both the inner and outer boundaries, by 
setting E = in the boundary cells. 

For our fiducial model we adopt para meters consistent with 
those observed for the RS Oph system (e.g., iHachisu & Katol2000l : 
ISchaef er 2009): Mi = 1.35Mo, q = M2/M1 = 0.6, R2 = 2.1 X 
lO'^cm, and P = 2n/n.2 = 453.5 days. (Note that these parameters 
are only self-consistent at the 10% level, due to observational un- 
certainties and rounding errors.) We assign a small initial surface 
density of E = Ig cm"^ at all radii, and allow the disc to be built-up 
self-consistently as material accretes from the secondary. 



2.2 S-curves 

It is well known that normal WD accretion discs, such a s those 
aroun d CVs, are thermally unstable at Tc ~ 10'*K (see, e.g.. lLasot3 
I200E and references therein), due to the rapid increase in opac- 
ity with temperature as hydrogen is ionized. However, it is also 
well-known that thermal instability alone cannot reproduce the be- 
haviour of CVs; instead we require the thermal-viscous instability, 
where the efficiency of angular momentum transport also increases 
substantially as hydrogen is ionized. This behaviour is incorporated 
in our model by assuming that a varies with temperature as follows: 




ffhot - Ccold 

l+exp(I^) 



+ ttcold • 



(16) 



100 1000 
I / g cm~^ 

Figure 1. S-curves for our disc model, computed for the parameters of 
RS Oph system. The upper panel shows the curves in the standard E-T^jf 
plane; the lower panel the more instructive E-M plane. In each panel the red, 
black and blue curves represent S=10"', 10"& lO'^cm respectively. The 
dotted and dashed curves are computed for fixed a (cold and hot, respec- 
tively), while the solid curve uses the temperature-dependent ^-prescription 
described in Equation 1161 (which "jumps" from ffcou on the lower stable 
branch to ffi,ot on the upper stable branch). 



This form allows for a smooth transition between the "cold-state" 
value cFcoM and the "hot-state" value o-hot, as the temperature T rises 
above some critical value Ta-it- In our fiducial model we set TctK = 
5000K and AT = 250K, with a^^ = 0.01 and tthot = 0.1. 

From these equations we can determine the range of param- 
eter space where the disc will be unstable by assuming local ther- 
modynamic equilbrium (Qy = Q-). At a given radius in the disc 
we can then solve for E as a function of T^, using Equations [811151 
above. As the discs are optically thick everywhere we can define 
the effective (or photospheric) temperature T^fs = TcT"'''', and de- 
rive the "S-curves" characteristic of disc instability models. These 
are shown in Fig.[T] 

The S-curves from our simple, one-zone model are some- 
what crude, and l ack the detail obtai nable with more sophisti- 
cated models (e.g.. lHameurv et ai]|l99 8). but they capture the es- 
sential ingredients of the model well. We immediately see from 
Fig.[T] that a sufficiently large disc is unstable over a very wide 
range in accretion rates: if the disc is ~ lO'^cm in size, then it 
will be unstable at some radius for all accretion rates in the range 
10"''MGyr"'< M < lO^^Moyr"'. Further tests indicate that the 
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shapes of the S-curves vary only sUghtly with changes the numeri- 
cal parameters T^f^^ or Ar, and that the general behaviour is robust. 



2.3 Emitted spectrum model 

In order to compare to observations, we compute the (continuum) 
emission from the binary system at each time-step. The model has 
three contributions: the emission from the disc, the emission from 
the boundary layer on the WD surface, and the emission from 
the secondary. The secondary is assumed to be a Roche-lobe fill- 
ing giant (with photospheric radius 0.7 times the Roche radius, 
R, - 120Ro), and is modelled as black-body with temperature T.. 
The emission from the disc is computed by assuming that each an- 
nulus emits as a black body with temperature T^g. To compute the 
emission from the primary we assume that half of the accretion lu- 
minosity is radiated in a boundary layer on the WD surface. If the 
boundary layer radiates as a black-body, we therefore have 



GMiM 



(17) 



where the accretion rate M is computed at the inner boundary. If the 
WD radius is known Equation[T7]uniquely specifies the bound- 
ary layer temperature Tbi- We adopt a radius R^d = 2.5 x lO^crn 
(consistent with a WD mass of 1.35Mo, e.g.. lAlthaus et alj|2005l) 
and a secondary surface temperature of T, = 3500K ( Schaefej 
With this parametrization, the total accretion luminosity 
scales as 



Lbi = 1.1 X 10" ergs"' 



Ml 



M 



1.35MgA5 X lO-'Moyr-i 

^/_^wd_r' 
\2.5 X lO'^cm/ 



(18) 



and the boundary layer temperature as 



7.1 X 10^ K 



Ml 



I.35Mf 



1/4 



M 



5 X lO-'Moyi- 



1/4 



^\2.5 X 10«cm/ 



-3/4 



(19) 



The luminosity of the secondary is L, = AttRIosbT* = 7.5 x 
10^* ergs"'. 

The inner edge of our computational grid is at = lO'^cm, 
and is thus an order of magnitude larger than the radius of the WD. 
It is very computationally expensive to extend the grid to smaller 
radii, because the innermost grid cells require very short time-steps 
during outburst. Neglecting this iimer part of the disc has a negligi- 
ble effect on the models, as the boundary plays no significant role 
in the evolution, but it can have a potentially important effect on 
the SED (as we neglect the hottest part of the disc). In outburst the 
disc follows an almost perfect T^g oc R^^/* power-law, so we also 
compute an additional SED in which we extrapolate the disc tem- 
perature beyond the inner edge of our computational domain. We 
extrapolate in to = 8 x 10**cm, using a T^g oc R^^^* power-law 
and normalising at R = 5 x lO'^cm. This extrapolation is some- 
what dubious, as departure from a simple power-law is likely in the 
boundary layer, but as it has only a small effect on the resulting 
light-curves (see Fig.|4]l we do not consider it to be a significant 
source of uncertainty. 




t / yr 

Figure 2. Accretion rates on to the WD as a function of time for the fiducial 
models (computed at the inner boundaiy of the grid). The black, red and 
blue curves represent the results for log[o(/Wi„fiiii/Moyr"')= -6.0, -7.0 & 
-8.0 respectively. The offsets between the outburst times for the different 
models is an artefact of our initial conditions, as the models with lower 
infall rates take longer to build up sufficient mass to trigger their initiai out- 
burst. Once the fimit cycie has been estabfished aif of the modefs show a 
recurrence time-scale of ^ 12yr Note aiso the pattern of aftemating fai'ge 
and smaff outbursts in the modefs with farger infaff rates. The hatched re- 
gion shows the approximate locus of the stabfe nucfear burning b and for a 
1.35MoWD, based on the cafcufations of lShen & BifdstenI j2007l) . 



3 RESULTS 

The model described above essentially has two free parameters: the 
rate of accretion from the secondary, and the angular momentum of 
the accreted material. These are parametrized in our models as the 
integrated accretion rate and the radial profile of the source term 
^infall- For simplicity, we assume that the profile is a (5-function at 
the circularisation radius R^., and thus the model is completely spec- 
ified by the two parameters MjnMi and R^. 

We assume that accretion from the secondary proceeds via 
Roche lobe overflow. In this case the circularisation radius can b e 
calculated explicitly from the binary parameters (e.g., Wvnn"2008!), 
and we find that R^ = 1.74 x lO'^cm. We compute the evolution 
of models for a range of infall rates log|o(Minfaii/Moyr"')= -8.0, 
-7.5, -7.0, -6.5, -6.0. In each case the disc builds up due to the ac- 
cretion of mass from the secondary, and once the disc is sufficiently 
massive we see repeated outbursts due to the thermal-viscous in- 
stability in the disc. The outbursts are of the "inside-out" type, and 
are triggered close to the inner boundary. After some initial tran- 
sients have died out a stable limit cycle is established, with the 
all of the discs showing repeated outbursts with a recurrence time- 
scale of ^ 12yr. Fig.|2]shows the accretion rates on to the WD as a 
function of time. The magnitude of the outbursts depends on how 
much mass has accumulated in the outer disc during quiescence, 
with higher values of Mi„faii resulting in much larger outbursts. The 
magnitude of any given outburst is determined by how far outward 
the ionization front can propagate, which in turn depends on the 
distribution of mass 'L(R) at the moment the outburst is triggered. 
For high values of Mi„faii we see a pattern of alternating large and 
small outbursts, with the larger outbursts propagating much further 
into the outer disc and accreting substantially more mass. Note, 
however, that in all cases the time-averaged accretion rate on to the 
WD is exactly equal to the rate of infall on to the disc from the 
secondary. 

The fact that the recurrence time-scale is independent of Mi„fau 
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R / cm 

Figure 3. Surface density evolution during an outburst in tlie model with 
log(()(Mi„faii/Moyr"')= -6.5, for both large (top) and small (bottom) out- 
bursts. The upper and lower solid curves respectively show the surface den- 
sity profiles immediately before and after the outbursts; the dotted curves 
show evolution during the outburst. The outbursts are triggered at small 
radii, and the ionization front then rapidly propagates outward thi'ough the 
disc. In the upper panel the curves are plotted 1.3, 1.7, 2.6 & 3.55yr after 
the initial profile; in the lower panel the corresponding times are 1.3, 1.5, 
1.95&2.45yr. 

can be readily understood. Fig.|3] shows the surface density evolu- 
tion of the disc during both large and small outbursts. Following an 
outburst the disc is essentially empty at small radii, but the structure 
at larger radii depends on how far the outburst was able to propagate 
into the disc. This in turn depends on the size of the mass reservoir 
that has built up in the outer disc (at radii ~ R^): if the outer disc is 
sufficiently massive then the ionization front can propagate to larger 
radii. We see from Fig.[T]that the disc is expected to be unstable at 
some radius for all of the infall rates considered here, so the recur- 
rence time-scale is simply the time required for the surface den- 
sity at some radius to exceed the critical value for instability. This 
first happens at small radii, typically a few times lO'^cm (hence the 
"inside-out" outbursts), so the recurrence time-scale is essentially 
the time required for gas to accrete from down to the triggering 
radius. This in turn is simply the viscous time-scale at R^, which 
in the cold state depends primarily on a^ou (with only a very weak 



Llghtcurves for fiducial models 
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Figure 4. Simulated V-band llghtcurves for the fiducial models (in absolute 
AB-magnitudes): note the small dynamic range on the vertical axis. As be- 
fore, the black, red and blue curves are for log[()(Mi„faii/M0yr"') = -6.0, 
-7.0 & -8.0 respectively. In each case the thinner line shows the Hght-curve 
if the disc is extrapolated beyond the inner edge of the grid (as described in 
the text); in practice, we see that this makes a negligible difference to the 
optical light-curve. The largest outbursts only marginally exceed the quies- 
cent flux level, and are barely detectable above the emission from the giant 
secondary. By contrast, during outburst s RS Oph is observed to brighten by 
6-7 mag in the V-band jSchaefeill2010l) . 

dependence on E, through the dependence of the viscous time-scale 
on the disc temperature). Larger infall rates therefore lead to larger 
and longer-duration outbursts, but do not alter the recurrence time- 
scale significantly. This is verified by computing additional models 
with ttcoid = 0.005: here, the recurrence time-scale is consistently 
=^ 25yr, almost exactly a factor of two longer than in the models 
with ci-coid = 0.01. 

Fig.|4]shows the simulated V-band light-curves for the fiducial 
models. Even in the largest outbursts the luminosity of the disc only 
barely exceeds that of the giant, and at optical wavelengths the net 
effect is a brightening only at the ~ 0.01 mag level. By contrast, in 
the V-band RS Oph is observed to bri ghten by a factor of 100-1000 
(6-7 mag) during its RN events (e . g . ~ I S chaef eilllO foH'l . The differ- 
ence between the "extrapolated" and standard light-curves from our 
models is essentially negligible, as the inner disc (the extrapolated 
part) contributes negligibly to the (/-band ffux. In the largest out- 
bursts the bolometric luminosity of the boundary layer (effectively 
the accretion luminosity) exceeds that of the secondary by almost 
two orders of magnitude, as can be seen from Equation [TS] How- 
ever, because of the high temperature of the boundary layer (Equa- 
tion|19t most of this luminosity is emitted at short wavelengths, and 
at optical wavelengths the outbursts are only marginally brighter 
than the secondary. 

During outbursts, the observed optical llghtcurves of RS Oph 
show a rapid rise (on a time-scale of days) followed by d ecay back 
to the quiescent level on a time-scale ~ 100 days (e.g.. ISchaefeil 
I2OIOI) . The total duration of the outbursts in our models is some- 
what longer than this, by a factor of a few (see Fig.[5j, but given 
the uncertainties in our model parameters (notably ofhot) is not clear 



The bolometric luminosity of the outbursts remains poorly constrained, 
primarily due to the large uncertaint y in the distance to RS Oph, but is es ti- 
mated to peak at ~ lO^^^ergs"' (e.g.. lSchaefe3l2009l ; IOsbome et alj201 ih . 
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Properties of outbursts 




Outburst duration / yr 

Figure 5. Duration and total accreted mass for 17 consecutive outbursts 
in each of the fiducial models. The smaller outbursts are comparable for 
all values of M|„f;,ii, but the larger outbursts increase sharply in magni- 
tude as Minfaii increases. The largest outbursts result in the accretion of 
> lO^^Moon to the WD on time-scales of 2-3yr The hatched region de- 
notes where the mean accretion rate during an outburst lies in the sta- 
ble nuclear burn i ng ba nd for a I.SSMqWD, based on the calculations of 
IShen & Bildsteij i2007l) . 



rence time-scale, and argues strongly against thermonuclear run- 
away as the origin of the observed outbursts. 

By contrast, Fig.[2] shows that in most cases the accretion rate 
during outbursts is sufficient to allow stable nuclear burning on the 
WD surface. For a 1.35MqWD the efficien cy of nuclear burn ing 
exceeds that of accretion by a factor of ~ 20 jprank et"al]|2002h . so 
comparison with EguationllSlshows that stable nuclear burning can 
yield a luminosit}|3 of ~ 10'*'ergs"'. This is comparable to the ob- 
served outburst luminosities, and shows that stable nuclear burning 
is an energetically plausible mechanism for powering the outbursts. 
However, in most outbursts the accretion rate exceeds the upper en- 
velope of the stable nuclear burning band, so it is not obvious that 
stable nuclear burning can be sustained. What happens in this case 
is far from clear. If the accretion rate exceeds the stable burning 
limit the accretion flow becomes thermally unstable and expands, 
in a manner similar to a red giant atmosphere, but the time-scale 
for this expansion depends sensitively on several unknown param- 
eters. It is possible that such an event c ould be responsible f or the 
observed outbursts of RS Oph (see also lKing & Pringl3l2009il . but 
detailed modelling of this process is beyond the scope of this paper. 
We note also that this type of mechanism, where disc instability acts 
as the trigger for the observed outbursts, may provide an explana- 
tion for the recently-detected pre -outburst signal in the lightcurve 
of RS Oph jAdamakis et"ai]|201 ih . 



whether this discrepancy is significant. However, the shape of the 
lightcurves in our model (which are essentially symmetric about 
the outburst peaks) is inconsistent with the observed lightcurves. 
Our model therefore fails to reproduce both the magnitude and time 
evolution of the outbursts that are observed at optical wavelengths, 
so disc instability alone cannot explain the outbursts seen in RS 
Oph. 

Fig.[5] shows the properties of the outbursts for the fiducial 
model set. For small values of Minfaii the outbursts are short in du- 
ration, and the WD accretes < 10 'Mq. However, the magnitude of 
the outbursts increases sharply with increasing Mmm^ and for val- 
ues Minfaii > 3 X lO^'Mgyr^'the total mass accreted during an out- 
bursts can exceed IO^'Mq. Also shown in Fig.|5](see also Fig.[2l( is 
the approximate locat ion of the stable nuclear burning band, based 
on the calculations of IShen & Bildstei] j2007h . Only for accretion 
rates within this region accreted hydrogen can burn stably on the 
WD surface: for lower accretion rates the WD simply accretes un- 
processed gas from the disc, while for higher accretion rates the 
WD surface becomes thermally unstable. We see that in our model 
the accretion rate on to the WD during outburst can span all three of 
these regimes, depending on the accretion rate from the secondary 
on to the disc. 

Given the similarity between the recurrence time-scale of the 
disc instability in our model and that of the observed outbursts, 
it is tempting to speculate that accretion due to the disc outbursts 
could be responsible for triggering a thermonuclear explosion on 
the disc surface. However, this requires the WD to accrete unpro- 
cessed hydrogen from the disc, while in most of our models the 
accretion rate exceeds that at which the accreted material will un- 
dergo nuclear burning as soon as it reaches the WD surface. In- 
deed, in all our models the mass of unprocessed material accreted 
by the WD (i.e., mass accreted below the nuclear burning band) is 
< IO^'Mq per outburst. This suggests that the time-scale required 
to accrete enough unprocessed hydrogen to trigger a thermonuclear 
runaway is > 200yr. This is much longer than the observed recur- 



4 DISCUSSION 

Our simple one-dimensional disc model has several limitations, the 
effects of which may be important in certain circumstances. The 
most significant of these is our almost total neglect of the disc's 
vertical structure, the so-called "one-zone" approximation (Equa- 
tion [T2). In general this is a reasonable approximation to make in 
thin discs, but it can break down in regions where the opacity varies 
significantly over the vertical extent of the disc. This can occur dur- 
ing disc instability outbursts (as the instability is triggered by the 
steep dependence of opacity with temperature), and an improved 
treat ment would consider the vertical structu re of the disc in detail 
(e.g.. Bell &Lirilll994 l Hameurv et al.ll993) . In practice, however, 
this procedure can be well approximated in one-dimensional mod- 
els by modifying the resulting S-curves appropriately, and the re- 
sults from such models are broadly consistent with our results. Thus 
while the precise details of our model thus depend on the shape of 
our derived S-curves [which in turn depend on the chosen functions 
Kip, T) and ji(p, T)], the qualitative behaviour should be robust. 

In addition our one-dimensional treatment explicitly neglects 
non-axisymmetric effects, in particular the tidal field of the binary 
and the localised nature of the accretion hot-spot. Our azimuthally- 
averaged approximation of the hot-spot heating (Eguationll It finds 
that the magnitude of the bulk heating is negligible, so while the 
hot-spot may have observational signatures the treatment of the hot- 
spot is unlikely to affect the disc evolution significantly. Further, 
because the disc instability outbursts are triggered at small radii, 
tidal effects (which primarily affect the disc near its outer edge) 
will have a negligible influence on the disc's evolution. 



^ Alternatively one can estimate this value as eMsteadyC^ - 1^3 X 
lO'^ergs"', where e = 0.007 is the efficiency of hydrogen burning. 
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4.1 Burning during outbursts? 

We have shown that the accretion disc in RS Oph is likely to un- 
dergo dwarf nova outbursts, with a characteristic recurrence time 
of order ^ 12 yr. There is a pattern of alternating large and small 
outbursts, with the large outbursts recurring at intervals of ^ 25yr, 
quite similar to the observed recurrence time-scale of RS Oph's out- 
bursts. We also find that even during the largest outbursts the pre- 
dicted accretion luminosity is insufficient to explain the observed 
optical outbursts of RS Oph. However, in all cases the accretion 
rate during these large outbursts crosses the band for steady nu- 
clear burning, and this offers a plausible explanation of the ob- 
served outbursts. Is it not clear how the system responds to both 
surface burning and a rapidly changing accretion rate, but our re- 
sult suggests that quasi-steady nuclear burning does occur at some 
point during these outbursts. Within the model used in this paper 
it is difficult to go further, as the development of the large out- 
bursts must be affected by the the high nuclear luminosity of the 
WD. The effect is probably similar to soft X-ray transients, where 
X-ray irradiation by the central accreting neutron star or black 
hole traps the disc in the hot state and prolongs the outburst un- 
til the central region s of the disc are largely drained of matter (e.g., 
iKing & R itter"l998^. The possibility that long-period dwarf novae 
outbursts m ight tri gger s t eady n uclear burning was considered by 
iKing, Rolfe, & Schenked ( I2003) . who showed that this may be a 
promising candidate for the progenitors of SNe la. If stable nu- 
clear burning can be sustained the WD mass will grow, and our 
model suggests that in this case the WD will be able to accrete 
> 10"'' Mq per outburst (see Fig. 5). 10''-10^ outbursts would there- 
fore sufficient for the WD mass to increase by 0.01-0. IMq, which 
suggests that if RS Oph is the progenitor of a Type la SN it could 
reach Mq, in < IMyr. 

This scenario, however, raises a number of questions about the 
long-term evolution of the system. If the WD mass is very close to 
Mch then a relatively small number of outbursts may be sufficient to 
trigger a SN. The WD mass has not been measured directly, how- 
ever, and if (as we argue) the RNe outbursts of RS Oph are not due 
to thermonuclear runaways there then the WD mass may be signif- 
icantly lower, ^ l.O-1.3M(70. In this case we require the outburst 
cycle to be maintained for ~lMyr in order to reach Mch, which in 
turn requires that the disc instability limit-cycle be maintained over 
a large number of outbursts. If a significant fraction of the WD 
mass is transferred from the secondary to the primary the binary 
will expand, and over time this may stop accretion via Roche lobe 
overflow. However, this so-called nova hibernation is usually short- 
lived, and typically evolution o f the secondary w ill re-start mass 
transfer relatively quickly (e.g.. lShara et al ] ll986h . The long-term 
stability of the outburst cycle is certainly an issue that must be con- 
sidered, but as long as the WD is significantly more massive than 
the secondary there is no obvious reason to expect binary evolution 
to prevent the WD accumulating mass on secular time-scales. 

In addition, repeated bursts of accretion in the stable nuclear 
burning band presumably result in a build-up of helium-rich ma- 
terial on the WD surf ace, and this could eve ntually lead to a he- 
lium nova explosion ( llben&Tutukovlll99lb . The accretion sce- 
nario considered here (burning on arrival) is rather different from 

^ It is also worth noting that if the system is to become a SN la then 
the initial (carbon-oxygen) WD mass m ust have been < 1.2Mo (e.g., 
Ilben & Tutuko\<1984lNomoto et al.ll984 . A larger initial mass results in a 
oxygen-neon-magnesium WD, which will collapse into a neutron star rather 
than exploding as a SN la (e.g. JSaio &Nomotoll"983) . 



that typically used in helium nova models, and it is not clear 
whether the conditions for helium ignition will be met before the 
WD is able to go supernova. However, detailed models of helium 
flashes on WDs find that helium novae eject signi ficantly less mass 
than they accumulate prior to the explosion (e.g.. lKato & Hachisul 
1 19991) . This is in contrast to conventional hydrogen novae, which 
typically eject all of the accreted mass, and suggests that the WD 
mass increases with time even if helium novae are triggered. Fur- 
ther study of this issue is still needed, but growth of the WD mass 
towards Mch through repeated disc instability outbursts seems to be 
a plausible path to Type la SNe. 



5 SUMMARY 

In this paper we have presented a one-dimensional time-dependent 
model for disc accretion in the RS Oph system. We find that the 
large disc around the white dwarf primary is always unstable to the 
thermal-viscous instability, and we are able to study the properties 
of these outbursts in detail. We find that: 

• Disc instability in the RS Oph system naturally leads to re- 
peated outbursts of accretion on ~ 10-20yr time-scales, similar to 
the recurrence time-scale of the observed outbursts. 

• The disc luminosity during outburst is not sufficient to explain 
the observed RNe outbursts. 

• The mass of unprocessed hydrogen accreted during the disc 
instability outbursts is small, < IO^'Mq per outburst. This argues 
strongly against thermonuclear runaway as the origin of the ob- 
served outbursts. 

• It seems possible that disc instability combined with surface 
nuclear burning can explain the observed RNe outbursts. The re- 
currence time-scale is naturally consistent with disc instability, and 
the luminosity of the observed outbursts is comparable to that of 
steady nuclear burning. 

• In our models the majority of outbursts attain (or exceed) the 
accretion rate required for steady nuclear burning. This suggests 
that the white dwarf is thus able to grow to the Chandrasekhar mass 
on <Myr time-scales, and we suggest that RS Oph and similar sys- 
tems may therefore be possible progenitors for Type la supemovae. 
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